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Stiitiiiiary 

A liiiite-dilFrrence scheme which uses fmiltipleoserset meshes to simulate the aerodynamics of aircraft, store interaction 
and store separation is <les<riheil. In this chimera, or multiple mesh, scheme, a complex conliguratinn is inap|ied using a 
major grid ahout the main conipoiient of the configuration, and minor overset meshes are used to map each additional 
component such as a store As a first step in modeling the aerod>namics of store separation, two-dimensional inviscid flow 
rah Illations have heen carried out in which one of the minor meshes is allowed to move with respect to the major grid 
Solutions of calibrated two-dimensional problems indicate that allowing one mesh to move w llh respect to another does not 
adversely affect the time accuracy of an unstead) solution. Steadv, inviscid three-dimensional computations demonstrate 
the rapahility to simulate complex configurations, including closely-pai ked multiple bodies. 

I. Introduction 

Today s htgh-s|ieed aircraft are intended to carry stores either externally or s^mi-submerged, and are rated as to 
their ability to deliver these stores to their targets reliably and accurately, without loss of aircraft perfornianre and agility. 
However, store-induced aerodynamic drag can signiPicantly downgrade aircraft performance, while aerody namic interaction 
mav cause the released stores to scatter, run into each other, or impart the aircraft It is necessary, then, to develop reliable 
methods to preilict the aeroflynaiiiirs of the store/ airframe interaction. 

It is difficult to test various carriage ronflgiiralioris and releases of stores in the windtiinnel, becau* stores cannot 
be freely released, «>r, if sting mounted, can not move with coinplele frertlom. Moreover, scaling effec‘< .rom the iiKwlel 
to true sixe are diflu iilt to overcome Flight testing is an option available for aircraft /store conngurations, hut without 
analytical aml/nr experimental results to relv on, flight teats can be unpredutable, expensive, and dangerous. 

the development of a iiutnerhal sitiitilaluui (apabilily using finite-difference leclinif|iies would aid in the design and 
testing of aircraft carrying store’s. A inafor obstacle to a coinputational approach, however, is the dillirultv in developing 
efferlive discretisation prcM esses for complex cotifigurahotts. While a single global mesh alKuil an ainraft with stores can 
lie generated, it is unlikely to be rornputatiooally eflicient, especially for viscous flow simulation. Moreover, allowing the 
stores to move would require a new grid to be generated after each time step. 

Ill lieu of using a single niesli for rniiiplex body configurations, multiple irii*slies ran be used to map the overall 
conrigiiration. The use of a multiple grid approach can yield belter grid resolution, simplify the applirali<iii of houndary 
roiiditiofis. and ea.se the task of grid generation. Some of the disadvantages to multiple meshes are that the flow solver 
must lie modified, and that the bookkeeping that tracks the relationships among the meshes ran he complicated. Two of 
the multiple mesh techniques currently emphasised in the literature are patched grids (refs. 1-3) and overset grids (refs. 
1«9). aithough both could be considered as special rases of a general multiple mesh approach. 

We have been working with a chimera grid scheme which uses ovcrf*t meshes. In this approach, minor grids are 
overset on a major grid to completely model a configuration. These minor meshes are generated somewhat independently 
and are overset on the major grid without requiring any special mesh boundaries. As such, the minor meshes ran l>e 
allowed to move with respect to the major grid. 

The purpose of this paper Is to demonstrate the feasibility of this chimera scheme for future application to the 
problem of airrraft/iiore Interaction and separation. Unsteady flow results have been rotnpuied in two fliniensUiiis, where 
one component in a multiple mesh system has been moved with respect to the fixed system. In three dimensions, steady, 
inviscid solutioni have been obtained for a generic wing/body /tail configuration, and for cloaely-packr<l multiple bodies 
that can ultimately be run as a TEH (triple ejection rack) configuration. 


A brief description of the overset chimera grid package is given in section tl. Section III contains information about the 
flow solvers used in both the two- and three-dimension applications. Details about the combined flow-solver /grid- package 
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codes are given in section IV. The results are shown in section V, and concluding remsrks follow in section VI. 
II. Chimera Grid Management Scheme 


Our mesh scheme is one In which a major grid stretches to the far field and is often generated abniii a main body 
element such as a wing. Minor grids are then overset on the major grid so as to resolve secondary features of (he 
configuration such as stores, nacelles, or Haps. In general, the minor grids are overset on top of the major grid withoni 
requiring the mesh boundaries to join in any special way. With the use of such an overset mesh system the task of grid 
generation is great Iv simplified because individual grids can be generated independently (provided they resoUe the How 
gradients) and then superimposed to form the overall mesh configuration. In this way, one can budd up a grid system to 
treat complex configurations and avoid severe mesh distortion. 

To illuslratc the chimera grid technique an airfoil with a detached Hap tucked underneath it is shcm ii rnrshed in figure 
1 using overset grids. Doth the airfoil and the flap are shown with body>conforming O-grids, with the minor Hap grid 
overset on the major, or airfoil, grid Since the flap is an impermeable body, the points of the major grid that fall within 
the flap (or a curve of the minor grid that circumscribes the flap) are excluded or blanked from the flow field solution 
The detailed view of figure I shows the poinia (marked with filled circles) on the major grid that are exr hided from the 
solution because of the presence of the flap. The boundary of this “hole" in the major grid l»ecoMies the first of two 
interface boundaries that arise because of the use of overset meshes. Data for this hole boundary must be supplied from 
the solution contained on the overset flap grid. The outer boundary of the minor grid forms th»» second interface in this 
problem, and nearby points from the major grid are used to supply flow values along the minor grid outer houndary. The 
flow firld solutions on each grid are connected by transferring boundary information across these two interfaces. 

In the current codes, each grid is independently advanced in time. The boundary information, however, is not 
transferred simultaneously, but is updated with each time step in the following manner. Vsing figure 2 for illustration, 
the flow field on the major grid is time advanced using boundary data at time level rt The minor grid solution is then 
advanced with the Rap surface and its outer boundary also held at time n. Then the flow variables at the hole boundary 
are updated with the information just solved for on the minor mesh (the dotted line on the minor grid). The effect of 
the Hap is now imposed on the airfoil solution. The flow variables at the outer boundary of the mirmr griff are likewise 
updated with the information from the major grid (the dotted line), so that the flap feels the effect of the airfitil. Kiiiativ, 
the body surface and far field boundaries are updated in the usual manner. 

All of the multiple grid bookkeeping (the determination of the interface boundaries and communication controls) is 
done automatically by the chimera grid management package. Once the grids are input to the code, the chimera grid 
package defines the relationships amoung the meshes, excludes or blanks out the points within another bod) to form holes, 
and finds the nearby points to l>e usml for the transfer of infomation across the interfaces. Details about the bookkeeping 
procedures can be found in references 6-9. 

tiecause of the overset grid technique, grids ran be generated about relatively complex foiifigiirations. nrcsiise of the 
automatic chimera grid inanangement package, one or more grids can move with respect to each other. 


III. The Flow Solvers 


The wersel mesh scheme has been combined with implicit approximate factorisation algorithms (refs. It), II. 12) ff>r 
solving the Kuler equations in both two and three dimensbns (refa. 6-9). However, ralcutations in whi^h one body moves 
with respect to another have only been carried out in tw<^dimensions. 


In generalised curvilinear coordinates and nondi mention alised variables, the three-dimensional Kuler equations in 
strong conservation law form can be written at 
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and th« metric terms etc are obtained from Xtit etc. ip the usual way. 

An implicit approximale-factorisation finite*difTerence scheme for the three dimensional rquations (I) can be written 
as 

(/ I h^^A^){I + hS^B^){l -f f 6„h 4 6^G ^) , (3) 

where A(j" - Q"** Q" and h = At for first-order time accuracy. The three-dimensional flux Jacobian metrics A , /) , C 
are deflned in reference 12. 


The spatial derivative operators \ id df are approximated with second-order acrurale cenlral finile-dilferriire 

operators. On the implicit side of equation ,3)« the second-order centra] differences produce block-tridiaKonal matrix 
operators (/ 4 kS^A^), (/ 4 h6f,n^), and (/ ♦ which must be inverted sequentially to obtain AC^" Additional 

details about the three-dimensional Euler code are given in references 9, 10, and 13. 

In this initial research phase of the developing of a two-dimensional unsteady moving body code, a diagonalixation 
version of the approximate factorisation algorithm was used in a deliberate trade off of time accuracy for reduced com- 
putational work. By diagonaliting the block tridisgonal operators into scalar operators, the resulting system can be more 
efficiently solved. The '‘diagonalised’* form of the algorithm is 

r<'|/ + w<Aj|Srj/ + A/i,A;ir,-'Ag" ^ «" , (D 

with 

Ae = A, = T; , (5) 

and y These matrices are shown in the references 10 and 13. This form of the algorithm is at best first 

order accurate in time, as the eigenvector matrices are functions of ( and tf. With the addition of numerical dissipation 
terms, this diagonalized algorithm reduces the block tridiagonal inversion to a series of 4r4 matrix rnulliplics and scalar 
tridiagonal or (depending on the dissipation terms, see reference M| pentadiagonal inversions. 

IV. Flow Solver (Tihiiera Interactions 


As mentioned in the introduction, a flow aolver must be modified to account for the use of multiple nieshea and the 
'^holes'* in the grids that result because grid points of one mesh can fall within the body boundaries of another. These 
hole points must be blanked or excluded from the flow field solution. 

Two minor modifications have been made to the flow solver to acommodate new boundary conditions and the blanked 
or excluded points in the meshrs. A tnodiilsr boundary condition routine was developed to set up the boumbir) valiH*s for 
each mesh, because different meshes in the same configuration might require different boundary treatment. Por example, 
the k I line in a major, or flwv field, grid might represent a free stream boundary, while the k I line on the minor grid 
might represent a solid body boundary. Interpolation and update routines were also added for the transfer of information 
across the various mesh interfaces . 


1 he main change in the flow algorithm itself is the treatment of the hole boundaries. The hole information from the 
chimera grid package h stored in an array, here called IBLANK, which is defined for each point on each grid as 


lUUNK ^ 


if a point is not hlankiHl; 
if a point is blanked. 


( 6 ) 


In the flow solver, each element in the left-hand-side matrix and right-hand-side vector hi multiplied b> the appropriate 
IIILANK value More the matrix is inverted, losing a aiinpte Sxft tridiagonal system for illustration with points I and r» 
blanked, the resulting equation set would look like 
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with no changes computed for AQ 4 or AQs. The blanked values and QJ are updated in the interpolation routine 
outside of the flow solver. Pentadiagonal mairket are treated In a similar manner. With this approach, no s|iecial routines 
or kigk testa are required to exclude the blanked points from the flow field solution. 


Unsteady Modifications. In steady-state appikalioiisof the overset grids, the flow solver and the chimera grid manage- 
ment package were run Independently. The grid package was run once, and the IBLANK array and transfer itiformallon 
arrays were passed on to the flow solver. For an unsteady, moving body case where one body k free to move relative to 
a second, the grid package must be called each lime the body and its grid move. For these unsteady rases, care must be 
taken to update the Interface boundaries in a time-accurate fashion. 
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In order to simuUte, for example, the flow about a wing with a detached moving store beneath it, the main routine of 
the unstendy code was constructed in a modular fashion with four main units that are called each time step These units 
are: the chimera grid management package, a flow solver routine adapted for multiple meshes, a trajectory routine (for the 
present, a simple, predetermined motion), and a package to update the interface boundaries caused by the overset meshes. 
Rach of these is run s^^mewhat independently, with interactions among them limited to passing arrays of informatitrn 
With this modular approach, it is easy to replace one unit with another, suitably modified program without affecting the 
rest of the unsteady code. 


Figure 3 illustrates the interactions among the four modules. At the beginning of each new time step, the trajectory 
routine is called to calculate the new spatial position of the minor grid The chimera grid management package redeflne^ 
the relationship of the grids by finding new blanked points and locating and storing new nearby points to update the 
interface boundaries. If necessary, new grid metrics are formed. The flow field solution for each mesh is then advanced to 
the new time level, n t I, and the interface boundaries are updated for the next lime step. To facilitate output of interim 
flow field solutions, all information at the start of the sequence is at the old time level, and, by the end of the sequence, 
all information has been updated to the new lime step values. 

V. RMulta 


The overset grid procedure has been used to compute two-dimensional inviscid flows about geometries where one 
body is moved with resprcl to the rest of the con figuration. It has also been applied to steady, inviscid three-dimensional 
complex configurations. To verify the two-dimensional overset mesh code in the time-dependent mode, the transonic flow 
about an airfoil undergoing an naciilating plunge was run. A multiple mesh case with a minor O-grid wrappe«l around a 
NACA 64A0I0 airfoil overset on a Cartesian major grid is compared to a. single mesh case with a larger ()-gri<l aboni the 
same airfoil. The airfoil and its mesh move with a velocity of 

y* - - Awco§(u;: ) , (fi) 


where 


A - 0.0436chord, k - .4, and uf ~ AHord . 


The maximum upward excursion of the airfoil from *ero is w 4?l of the airfoil chord. Time histories of the lift coeflicient 
are shown in figure 4. Each run was started from steady slate (Mo© = O.fi, a ~ O"). The 720 step per cycle single mesh 
curve shown in figure i is considered to be the correct range for (ref. 15). A comparison between the single and 
multiple mesh time histories shows that the two are in phase, but that the multiple mesh results have a slight error in 
amplitude. This error in lime amirary can be improved by decreasing the time step, and is probably due to lime lagging 
the outer boundary of the minor grid and possible errors in interpolation. The error is small, and the results demonstrate 
the capability of the unsteady code to simulate the flow over a configuration with one component moving. 


Two transonic store separation simulatkins using predeterinined paths were computed with the Iwo-dimetisiotial 
unsteady code. Ilolh computations were started from the same steady-state solution. The initial grid arrangement for the 
store separation sinmlattons is shown in figure 5. Both bodies are meshed with body -conforming O-grids. The airfoil body 
is a thick ellipse with a chord length of I. The minor store body is a ellipse reduced by Filled circles mark 
the blanked {mints forming the hole in the major grid due to the minor bcwly and the hole in the minor grid caused b\ 
its intrrsi'ction with the airfoil. The empty circles denote the nearest points on the major grid usetl to update the minor 
grid outer boundary. The minor iMwfy ntrse is initially positioned at e/e - .3 and y/c * .15. The How field solution 

for the initial position is illustrated with Mach contours and is shown in figure 6. The free stream Mach number for this 
case is M™, ' -7, and lioth bodies are at 0” angle of attack. In the common overlap region, within the outer boundary of 
the minor grid and the hole in the major grid, both sets of contours are shown. The contours at the hole boumlaries and 
the outer boundary of the minor grid match well, showing that the communication between the grids across the interface 
boundaries is quite good. The flow channeled between the two bodies behaves somewhat like flow in a Iwo-dttnrnsional 
notsle, and there is a strong shock between the major and minor bodies, as well as some spillover on the lower surface of 
the leading edge of the store. There is some mismatch of contours at the shock caused by the interpolation between the 
two dissimilar meshes at the hole boundary. The slight discrepancy between the contours in the common sreas of the flow 
field is pariaitty attributed to the plotting package, as again the contoura are being calculated on two dissimilar meshes. 


The first store simulation case has a predetermined downward translation and rotation. The downward rale of .01 
corresponds to a drop of .75 rhordlengihs per second for an average realtime store separation profile. Ttir {>ilchdown 
rotation on the store is appronimately 2" per 100 time steps (A<Ooo/c - 0.025). Figures 7 and 5 show the grid positions 
and Mach contours for the store case after 400 time steps. The store body has dropped about 10^ of chord and has an 
angle of attack of fl^. Again, the filled circles mark the blanked points and the open circles denote the nearest points 
for update of the minor grid outer boundary. In figure S, the Mach contours reflect the movement of the store. A strong 
shock occurs on the lower surface of the store due to the downward pilch to the freestream, and the shock between 
the two bodies Is mtKh weaker. The mismatch of contours at the outer boundary of the miner mesh Is attributed to the 
interpolation and the lime lag on that boundary. Figures 9 and 10 show the grid positions and Mach contours for this 
cate after 600 lime steps. The store haa dropped another of chord and Is now at a 12* angle to the freestream. The 
shock between the two bodies has deteriorate to a weak shock at the tall on the upper surface of the store; it no longer 
spans the distance between the bodtea. 
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The second store separation simulation involves a predetermined path in which the store first drops down, then rotates 
hark up. if the caiculation were continued, the store would hit the wing. The first two figures, 11, and 12, show the grid 
postions and Mach contours after 100 time steps (with the same time step as before: AfOno/** ' 0.02J>) The store is 
dropping at the same rate as for the first case, with a downward rotation of the trailing edge of T per 100 time steps. 
The flow field (figure 12) shows the strong shock between the two bodies. After the 400th time step, the store starts 
to move upward at the same rate as before, but the trailing edge of the store continuing to rotate downward. The grid 
positions after 600 time steps are shown in figure 13, with the corresponding flow field solution shown in figure 11 The 
shock between the bodies is stronger than before and has moved forward. 


Steady, inviscid three-dimensional flow about a generic wing/body /tail configuration has been computet) and compared 
to experimental data The configuration is comprised of an ogive-cylinder fuselage with a wing and tail formed by swept, 
untapered wings with NAC'A 0012 airfoil cross-sections. Four grids are used to mesh this configuration Figure ir> shows 
the configuration with the surface grids for the fuselage, wing and tail, along with the outer boundaries for the wing and 
tail meshes. These outer boundaries did not intersect. The outermost of the four grids is used to resolve the solution from 
the near- to the farfield and is a warped hemispherical shell whose polar axis is coincident with the fuselage centerline. 
The me^h contains 37,000 points (74 x 2h x 20) and extends from 9 to $1 radii from the fuselage. The fuselage mesh is the 
second mesh in this configuration, another warped, hemispherical shell whose inner boundary is the fuselage surface. The 
grid contains 77,700 points distributed as 74 x 35 x 30 and extends to 11.5 radii. The outer boundary of the fuselage grid 
overlaps the inner boundary of the outer mesh by 2.5 radii. The wing and tail grids are warped, cylindrical meshes whose 
axes are directed along the wing- and tailspans. The end surfaces containing the roots of the wing and tail are coincident 
with the fusriage surface. The wing mesh has 27,720 points (66 x 211 x 15) with 20 spanwise points on the wing surface, 
and the (ail grid contains 15,120 points distributed as 56 x 18 x 15 with JO points along the tailspan. The composite mesh 
consists of 157.540 points. A solution was computed for this configuration at a Mach number of 0.9 at an angle of attack 
of 2*^. Figure 16 shows the comparison of experimental data (ref. 16) and the computation. 


As an example of another three-dimensional steady-state simulation, the flow about a generic three- body TKK-like 
conflguration was computed. The configuration consists of three elli|Mnidal bodies in a triangular arrangement (figure 17) 
The grids of the two small bftdies have major and minor axes one half those of the larger body. The two smaller bodies are 
embedded in the mesh of the larger, as indicated in figure 17. The figure slums the surface grids for all thri'e bodies, and 
the warped outer shells of the two minor meshes. All the grids are spherical. The mesh of the large ellipsoid has 26.250 
points (30 X 35 x 25); the two smaller meshes each have 15,750 points (,10 x 35 x 15). The composite mesh has 57,750 
points. Surface Mach contours, computed at = 0.8 and o = -2**, are shown in figure 18. 

VI. Concluding Reiniirka 

The chimera overset grid technique has been suiressfully demonstrated for the simulation of flow about complex con- 
figurations. While restricted here to two dimensions, the feasibility of numerically simulating transonic store srparation has 
l»een carried out. Moreover, preliminary steady-stale simulations about generic wing/body /tail and TKIf configurations 
show that the technique can be applied in three dimensions as well. At least for simple configurations such as the TFH, 
the simulation of unsteady three-dimensional store separation appears feasible in the near term 
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Figure 1 . Overset grids for an airfoil wi'h flap. 
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